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Introduction 

This report details work carried out under NASA contract NAS8-38609 delivery order no. 79. 

, A brief description of enhancements made to the NASA MSFC coherent lidar model is provided. 
Notable improvements are the addition of routines to automatically determine the 3 dB misalign- 
ment loss angle and the backscatter value at which the probability of a good estimate (for a maxi- 
mum liklihood estimator) falls to 50%. The ability to automatically generate energy/aperture 
parameterisation (EAP) plots which include the effects of angular misalignment has been added. 
These EAP plots make it very easy to see that for any practical system where there is some degree 
of misalignment then there is an optimum telescope diameter for which the laser pulse energy 
required to achieve a particular sensitivity is minimised. Increasing the telescope diameter above 
this diameter will result in a reduction of sensitivity. These parameterisations also clearly show 
that the alignment tolerances at shorter wavelengths are much stricter than those at longer wave- 
lengths. 

A brief outline of the NASA MSFC AEOLUS program is given and a summary of the lidar 
designs considered during the program is presented. A discussion of some of the design trades is 
performed both in the text and in a conference publication attached as an Appendix. 
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1.0 Improvements to the Lidar Model 

The following improvements have been made to the NASA MSFC coherent lidar model [1]. 

1.1 Back solution of misalignment angle required for 3dB loss 

In order to assess the alignment tolerance of different lidar point designs, it is convenient to look 
at the misalignment angle required to give a 3 dB SNR loss. Code to automatically determine the 
misalignment angle required to achieve a 3 dB misalignment loss was added to the existing model 
and a button added to the tool bar to enable easy user access to this function. This means that the 
user can now change the telescope diameter and then have the program calculate the correspond- 
ing 3 dB misalignment angle from any where within the model. The angle is calculated by back 
solving from the required solution to determine the input variable. The code is listed in Appendix 
(la). 


1.2 Back solution of backscatter required for 50% probability 

At low SNR values the velocity estimation algorithms may incorrectly identify a noise spike as 
the signal returned from the atmosphere. The probability that such an event will occur is directly 
related to the sensitivity of the instrument. The coherent lidar community has adopted the value of 
backscatter at which this probability is 50% (i.e. 50% of the time the estimator correctly estimates 
the wind velocity) as a measure of the sensitivity of a coherent lidar system and all the AEOLUS 
designs discussed in this report use this as a sensitivity comparison point. Code was developed so 
that when a new instrument design has been entered into the model a toolbar button can be used to 
initiate an automated back-solution to find the backscatter value required to achieve this sensitiv- 
ity value. The code is listed in Appendix (lb). 

1.3 Energy-aperture design trades code 

The energy-aperture product (EAP) has traditionally been used as a figure of merit for a lidar 
design and plots of pulse energy vs. telescope diameter for constant atmospheric backscatter have 
been used to determine design points for coherent lidars [2]. Such plots assume perfect alignment 
of the optical system and for a fixed pulse energy show continuous improvement in sensitivity as 
the telescope diameter is increased. Unfortunately in any practical system, particularly those 
involving rotating optical components such as the scanning system employed in coherent lidars, 
there is the potential for an angular misalignment to occur between the received beam (signal) and 
the local oscillator (lo) beam. Any angular misalignment in the system between the signal and lo 
beams results in a loss of mixing efficiency [3] and reduces the sensitivity of the instrument. The 
lidar model was enhanced by developing code (Appendix (Ic)) to plot the EAP as a function of 
misalignment angle for a fixed sensitivity (backscatter) or as a function of varying sensitivity with 
a fixed misalignment angle . Figure (1.1) shows EAP plots for both 2 pm and 9 pm lidars in which 
the misalignment angle has been allowed to vary whilst the 50% probability (Section 1 .2) for each 
wavelength was held constant at the value indicated. Figure (1.2) shows a similar plot for both the 
2 pm and 9 pm wavelengths when the misalignment angle is held constant at the value required 
for a 3 dB loss with a 0.5 m diameter telescope and the sensitivity is allowed to vary.Table (1.1) 
lists the model parameters used to generate the plots. 
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Figure (1.1) EAP plots for a 2 pm (top) and a 9 pm (bottom) lidar as a function of 
misalignment angle for a fixed sensitivity (see text). 

It can be seen that the misalignment sensitivity of the 2 pm system is much greater than for the 
9 pm system. It is easy to understand this by considering that the heterodyne mixing is dependent 
on the signal and local oscillator overlap function [3]. As the angular misalignment between the 
local oscillator and signal beam increases the degree of overlap also decreases. For a given tele- 
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Parameter 

Value 

Parameter 

Value 

Orbit height 

350 km 

Receiver type 

Complex 

Orbit inclination 

98 deg. 

Receiver geometry 

Wang 

Laser wavelength 

2.065479 pm 
9.1145187 pm 

Detector quantum effi- 
ciency 

0.6 for 2 pm 
0.4 for 9 pm 

Laser pulse energy 

See plots 

Scan mechanism 

Wedge 

' Laser pulse length 

0.5 ps 

System margin 

0.5 

Laser P.R.F. 

10 Hz 

Atmospheric model 

Midlatitude 
summer clear 

Telescope diameter 

see plots 

Target aerosol altitude 

1 km 

Nadir angle 

30 deg. 

Processing bandwidth 

± 20 m/s 

Transmit optics transmis- 
sion 

0.9 

Receive optics transmission 

0.9 

Polarisation efficiency 

1 

Misalignment angle 

see plots 

Tabl 

e (1.1) Parameters used for EAP plots. 


scope size the diameter of a diffraction limited beam at 2 pm, after propagation through a large 
distance, will be smaller than that for a 9 pm beam. Thus the overlap between the signal and local 
oscillator beams will decrease more rapidly at the shorter wavelength. 


It can also be seen that for a given misalignment angle there is an optimum telescope diameter at 
which the laser pulse energy required to achieve a particular sensitivity is minimised. This is very 
important for the design of a space-based lidar as the delay between output pulse transmission and 
signal return can be several milliseconds. If there is jitter in the spacecraft pointing during this 
round trip period it will be equivalent to an angular misalignment of the optical system such that 
even a perfectly aligned optical system will be unable to reap the benefits of a large telescope 
diameter. Thus thd spacecraft jitter, rather than tolerances in the optical design, may be the domi- 
nant mechanism limiting the choice of telescope diameter. Given the small misalignment 
(3.64 prad) required for a 3 dB loss at 2 pm for a 0.5 m telescope and the generally ill— defined 
stability specifications of many small spacecraft any design placing a short wavelength coherent 
lidar on a small spacecraft must be fully aware of this potential problem. 

1.4 Other Improvements 

A number of small improvements were made to the model. These included expanding the lookup 
table for the maximum liklihood estimator to cover a wider range of values of Omega and M, 
changing the misalignment loss vs. angle plot to an alignment efficiency vs. angle plot (more intu- 
itive this way), moving the “update performance” plots button to the toolbar so that it is accessible 
from everywhere, adding a few additional calculated fields and several other minor changes. 
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Figure (1.2) EAP plots for a 2 pm (top) and a 9 pm (bottom) lidar as a function of backscatter 

for a fixed misalignment angle (see text). 


5 







1.5 References 

[1] Gary D. Spiers, “Lidar Performance Analysis”, Final Report on NASA Contract No. NAS8- 
38609, Delivery Order 57. 

[2] See for example Figure 3-21 on page 3-18 of Lockheed Missiles & Space Company, “Phase II 
Design Definition of the Laser Atmospheric Wind Sounder (LAWS)”, DR-20 Vol. II: Final 
Report, Nov. 1992, Prepared for NASA George C. Marshall Space Flight Center, AL 35812 under 
NASA contract NAS8-37590. and also page 198 of TRW Applied Technology Division, “LAWS 
Laser Atmospheric Wind Sounder Two-Micron LAWS System Study”, Final Report, October 
1993, under contract NAS 1-19291 MSOA, Task 13. 

[3] Rod G. Frehlich, “Heterodyne efficiency for a coherent laser radar with diffuse or aerosol tar- 
gets”, J. Mod. Optics, 41, 11, 2115-2129 (1994). 
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2.0 AEOLUS 


AEOLUS is an acronym for Autonomous Earth Observing Lidar Utility Sensor and is a series of 
studies carried out at the NASA Marshall Space Flight Center towards developing a small space- 
based coherent lidar. A fairly complete description of AEOLUS and some of the design issues 
involved were presented in the paper “Direct global measurements of tropospheric winds employ- 
ing a simplified coherent laser radar using fully scalable technology and technique” presented at 
SPEE’s International Symposium on Optical Engineering in Aerospace Sensing, Technical Con- 
ference 2214 on Space Instrumentation and Dual-Use Technologies, session on “Faster, Cheaper, 
Smaller Space Science Optical Instruments,” Orlando, FL (6 April 1994). This paper is included 
here as Appendix II and is the main source of information concerning the AEOLUS designs. 

2.1 Description of the AEOLUS Designs 

The following tables summarises all of the AEOLUS designs analysed. 



Scenario 


Parameter 

1 

2 

3 

4 

5 

Units 

Wavelength, X 

9.11 

9.11 

2.065 

2.065 

2.065 

pm 

Pulse energy, E 

400 

400 

200 

200 

25 

mJ 

Pulse length, X 

0.5 

0.5 

0.2 

0.2 

0.5 

ps 

Max. pulse rate, PRF 

20 

20 

10 

10 

50 

Hz 

Detector Q.E. 

0.4 

0.4 

0.6 

0.6 

0.6 


Telescope diameter, D 

0.5 

0.5 

0.5 

0.5 

0.25 

m 

Scanning method 

Line scan 

Wedge 

Line scan 

Wedge 

Wedge 




Scenario 


Parameter 

6 

7 

8 



Units 

Wavelength, X 

2.065 

2.065 

2.065 



pm 

Pulse energy, E 

25 

25 

200 



mJ 

Pulse length, X 

0.5 

0.5 

0.2 



ps 

Max. pulse rate, PRF 

15 

15 

10 



Hz 


0.6 

0.6 

0.6 




Telescope diameter, D 

0.25 

0.233 

0.5 



m 

Scanning method 

Wedge 

Wedge 

Wedge 





Table (2.1) The AEOLUS configurations considered. 


Designs 1-4 were the initial point designs in which both 9 pm and 2 pm lidar systems were con- 
sidered. It was felt that the large rotating telescope used in the original LAWS concepts was too 
unwieldy for a small satellite design and so two alternate scan geometries were developed. The 
first used a pair of telescopes to point fore and aft to one side of the satellite (designs 1 & 3) - this 
produces a single line scan sampling pattem[l]. The two telescopes share a common volume and 
a common secondary mirror. The second concept uses a refractive wedge mounted in front of a 
nadir pointing telescope to produce a conical scan (designs 2 & 4) [ 1 ] . 

The dual-telescope approach was discarded because of limited coverage and awkward packaging 
and the 2 pm wedge scanned lidar (design 4) was favoured over the 9 pm equivalent (design 2) 
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because the need to cool the 9 pm detector introduced an additional power drain and thermal heat 
source (mechanical cooler) which has the potential to introduce vibration into the optical system. 
The 2 pm design also held the potential for the laser to be smaller and more compact than the 
9 pm equivalent. It should be noted that the pulse energies for the 2 pm and 9 pm designs listed in 
Table (2.1) were based on likely availability of a suitable laser from a vendor rather than from any 
lidar sensitivity issues and so they should not be assumed to be of equivalent performance. 

Design 5 attempted to develop a lidar small enough to fit on a Pegasus class launch vehicle. The 
limited sensitivity of such a small instrument indicated that most of the useful signal returns 
would probably be obtained from clouds. Therefore two backscatter values were determined for 
this instrument, one for a 1 km thick diffuse target (thin clouds or high aerosol concentration 
boundary layer) and one for a pulse width limited target (thick clouds). The instrument was 
designed to operate at a high PRF (50 Hz) to give a good shot pattern density and provide good 
statistical sampling of the atmosphere. Spacecraft accomodations at NASA MSFC determined 
that a combination of power and weight limitations aboard likely spacecraft limited operation of 
the instrument to ~ 30 % of each orbit. Discussions with the former LAWS Science Team indi- 
cated that they would prefer a greater orbit coverage and lower shot density. This led to design 6 
which is identical to design 5 with the exception of the PRF and orbit duty cycle, however the 
total number of shots per orbit remained the same. 

One consequence of operating the lidar over a whole orbit (design 6) was that ancillary equipment 
such as the receiver electronics also had to remain powered on. When this was factored in to the 
power budget the spacecraft accomodations people at NASA MSFC felt that the power budget 
was too high and so a design trade was carried out to arrive at design 7. 

Finally design 4 was revisited and the orbit duty cycle increased to 100 % in accordance with the 
preferences of the former LAWS Science Team. This resulted in two designs, one for a smallsat 
(design 7), Pegasus class launch and one (design 8) for a larger (but still small c.f. LAWS) instru- 
ment. 

2.2 Instrument Performance 

The basic parameter used to characterize the sensitivity of each design was the value of backscat- 
ter at which the probability of correct detection of the return signal from the atmosphere was 50%. 
This was calculated using the lidar model outlined previously[l]& Section 1.0. Table (2.2) and 
Table (2.3) list the common parameters used for the analysis. .whilst Table (2.4) summarises the 


Orbit height 

350 

km 

Receiver type 

Complex 


Orbit inclination 

98 

deg 

Receiver geometry 

Wang 


Nadir angle 

30 

deg 

Target altitude 

1000 

m 

Transmit optics efficiency 

0.9 


Horizontal wind velocity processing 
search space 

±20 

m/s 

Receive optics efficiency 

0.9 

H 

Vertical range integration (Clear Air) 

1000 

m 



■ 

(Clouds) 

63.67 

m 

Polarisation efficiency 

1 


Mid-latitude summer -clear atmosphere 


Table (2.2) Parameters common to all scenarios. 
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System margin 

0.5 


WPL-37 refractive index structure constant profile 


Table (2.2) Parameters common to all scenarios. 



Scenario 


1 

2 

3 

4 & 8 

5 & 6 

7 

Units 

Effective telescope diameter, D eff 

0.5 

0.433 

0.5 

wamm 

0.217 

0.202 

m 

BPLO/signal misalignment angle, a 

13.91 

16.07 

3.15 

3.64 

7.28 

7.81 

|irad 


Table (2.3) Effective telescope diameter and misalignment angle tolerance for each of the 
five scenarios assuming a nadir angle, 6 n of 30 degrees. 



Design No. 

1 

2 

3 

4 & 8 

5 & 6 

5b & 6b 

(3(50%) (/(m-sr)) 

@ 1 km vert, resoln. 
@ 64 m vert, resoln. 

7.7 x 10~ 9 

1.0 X 10" 8 

4.8 x lO " 8 

6.4 x 10- 8 

2.1 x 10~ 6 
2.3 x 10~ 5 

2.4 x 10~ 5 
2.6 x 10~ 5 


Table (2.4) AEOLUS backscatter values for a 50% probability of correct detection. 

sensitivity obtained from each system. 

2.3 Mass and Power Trade on Design 6. 

As was mentioned earlier, the AEOLUS design no.6 was found to be too power hungry (for a 
Pegasus class launch) by the spacecraft accommodations people. They provided an estimate that 
for every Kg of mass given up by the instrument, approximately 20 W of power could be gained 
from the spacecraft by increasing the solar array size. 

For the initial instrument design, the spacecraft could supply 200 W of power, which was insuffi- 
cient for the instrument, however it was realised that by reducing the diameter of the Ge wedge 
scanner, the weight of the instrument could be reduced thereby freeing up mass for the power sys- 
tem. Figure (2.1) shows the wedge mass as a function of its diameter. It can be seen that large 
weight savings can be obtained for small changes in the wedge diameter. Figure (2.2) shows the 
laser and instrument power as a function of laser pulse repetition frequency (PRF). Also plotted is 
the wedge optical diameter as a function of the laser PRF. 

For laser PRFs below 5 Hz the instrument power is less than 200 W and no weight reduction is 
needed, however above 5 Hz the power required is greater than 200 W and the instrument mass 
needs to be reduced to increase the size of the spacecraft power system. Thus above 5 Hz PRF the 
optical diameter of the wedge reduces. It has been suggested (Dr. G.D. Emmitt - private commu- 
nication) that for optimum sampling of the atmosphere with this instrument, a laser PRF of -10 - 
20 Hz would be advantageous. It can be seen that at a PRF of 15 Hz, the instrument requires 
-230 W. This power can be achieved by reducing the wedge optical diameter to 24.3 cm from the 
original 25 cm diameter. This results in a small 0.7 dB loss in SNR. Thus by trading weight for 
power the instrument design has been improved. It should be noted that this is a conservative 
trade, as neither the telescope primary mirror nor the support structure have been reduced . 
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Power (W) 



Wedge mass as a function of diameter 



20 21 22 23 24 25 26 

Physical Diameter (cm) 

(Optical diameter is 1cm smaller.) 


Figure (2.1) Dependence of wedge mass on diameter 


Trading wedge diameter for power 



PRF (Hz) 


Figure (2.2) Trading wedge diameter for increased power to the instrument. 
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Appendix I Lidar Model Code 


The code listed here follows the same conventions as listed in the previous report[l]. The code 
makes use of the OPTIMIZER routine built into Quattro Pro and it has been found that this rou- 
tine is very reliable when solving for a value if the initial guess is less than the final value. When 
the guess value is larger than the final value the routine will occasionally fail to find a solution - 
thus each use of the OPTIMIZER explicitly sets the initial guess value for the unknown to a value 
much smaller than any solution likely to be encountered. 

Appendix (la) Code to calculate 3dB misalignment angle 


The following table is a coded macro to backsolve for the misalignment angle required to give a 
3 dB SNR loss. It is executed when the user clicks on the “3dB Loss” button on the toolbar. 



A 

B 

c 

a 

THIS MACRO CALCULATES THE BPLO/SIGNAL MISALIGNMENT ANGLE REQUIRED FOR A 3dB 
LOSS 

38 




39 




31 




a 

3dB_angle 

18.20474 


a 


7.81320 


a 


14.37397 


a 




a 




a 




a 

3dB_calc 

(LET $Current_angle, @ROUND($Align_angle,6)} 

Store current angle 

a 


(LET $Current_loss_dB, @ROUND(- 
10*@LOG($eta_misalign),6)} 

Store current misalignment angle 

m 

m 


{Optimizer.Reset}{Optimizer.Solution_ceH 
[lidar]$eta_misalign } 

Set OPTIMIZER to solve for mis- 
alignment efficiency 

m 

■ 

- 

{Optimizer.Solution_goal Target Value} 

OPTIMIZER to aim for a target 
value 



{ Optimizer. Target_Value 0.5 } 

Set target value 

a 


{ Optimizer. Variable_cells [lidar]$Align_angle } 

Vary misalignment angle 

53 


{ Optimizer. Add 1 Align_angle" ,>=,"0" } 

Keep misalignment angle >=0 

54 


{Optimizer.Solve} 

Solve! 

56 





Appendix (lb) Code to calculate P(50%) 

This table contains macro code which is used to backsolve for the backscatter value at which the 
probability of a good estimate is equal to 50%. It is executed when the user clicks on the “P(50%) 
button on the toolbar.” 


n 























A 

B 

c 

SM 


mm 

U 

THIS MACRO CALCULATES THE BACKSCATTER REQUIRED TO GIVE A 50% PROBABILITY 

a 

OF A GOOD ESTIMATE AND INSERTS THIS BACKSCATTER VALUE INTO THE MODEL. 
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91 




a 

50%BETA_CALC 

{LET $Save_beta, $beta} 

Save existing value of backscatter 
[the save_beta variable is with the 
performance calculation macro) 

mm 

■ 


{LET $beta,le-13} 

Set the initial guess for beta to a 
very small value. 

RH 


{Optimizer.Reset}{Optimizer.Solution_cell 
[lidar]$Pgood } 

Set OPTIMIZER to solve for 
Pgood. 

■ 


{Optimizer.Solution_goal Target Value} 

Set OPTIMIZER to solve for a tar- 
get value of Pgood. 

91 


{ Optimizer. Target_Value 0.5 } 

Set target value of Pgood as 0.5. 

mm 


{ Optimizer. Variable_cells [lidar]$Beta} 

Tell OPTIMIZER to vary backscat- 
ter to achieve target value. 

91 


{Optimizer.Solve} 

Solve! 


Appendix (Ic) Code to calculate and display an EAP plot 

This macro is executed when the user clicks on the EAP button on the toolbar. When run it moves 
to a screen displaying a photograph of hurricane Gladys. This is purely for aesthetic reasons and 
provides the user with something to look at while the code is run as this can take some time if a 
large parameter space is selected. A dialog box is then displayed giving the user a choice of two 
types of EAP plot as discussed previously(Section 1.3). After a selection has been made another 
dialog box is then displayed to allow the user to enter the minimum and maximum values required 
for each of the independent variables and the number of intermediate steps between these 
extrema. The telescope diameter is treated as an independent variable for both types of EAP plot 
whilst the misalignment angle and backscatter alternate as the second independent variable 
according to the plot type selected by the user. When backscatter is an independent variable the 
intermediate values are calculated for a log axis whilst a linear scale is used for intermediate val- 
ues of misalignment angle. The value of the fixed parameter is taken from the current value in the 
model and the routine then solves for the pulse energy required to achieve Pgood = 50% for each 
combination of telescope diameter, backscatter and misalignment angle. When all the values are 
computed the routine edits the EAP plot, changing the series and updating axis labels and scalings 
as necessary to provide a reasonable plot. .The plot is then displayed to the user. 
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B 


C 



THIS MACRO CALCULATES AND DISPLAYS AN EAP PLOT AND ALLOWS THE USER TO 
CHOOSE 


77 BETWEEN CONTOURS OF CONSTANT BACKSCATTER OR MISALIGNMENT ANGLE. 


Rows 78 - 84 are all blank and have been skipped in this table. 


85 

eap_option 

86 

save_angle 

87 

save_tscope 

88 

save_energy 

89 

tscope_min 

90 

tscope_max 

91 

tscope_nstep 

92 

eapvar_min 

93 

eapvar_max 

94 

eapvar_nstep 

95 

loopl_count 

96 

loop2_count 

97 

row_number 

98 

col_number 

99 

tscope_step 

100 

eapvar_step 



16.0667307422832 


0.5 


Used to determine type of EAP plot 
required: 

1 = vary misalignment angle 
0 = vary backscatter 




eapvar is either misalignment angle or 
backscatter 



($TSCOPE_MAX-$TSCOPE_MIN)/ 

$TSCOPE_NSTEP 


($EAPVAR_MAX-$EAPVAR_MIN)/ 

$EAPVAR_NSTEP 


Counter for tscope diameter 


Counter for angle or backscatter 


Used to locate position in results table. 


Used to locate position in results table. 


Calculates incremental step to evenly 
space telescope diameter values on x-axis. 


Calculates incremental step for misalign- 
ment angle only NOT backscatter. 
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B 


C 



103 marker_style 


104 guess 


105 old_nstep 


@LOG($EAPVAR_MIN) 


( @ LOG($E AP VAR_M AX)-$ A)/ 
($EAPVAR_NSTEP) 


@IF(LOOPl_COUNT>9,+"M"&@STRI 
NG($LOOP 1 _COUNT,0) &" ,5 ,0" ,+"M0" 
&@STRING($LOOP1_COUNT,0)&",5, 
0 ") 


a and bval are used to uniformly distribute 
backscatter values alnog a logarithmic 
scale. 




This expression constructs a marker style 
tag so that each data series has a unique 
symbol. 


Used to set up OPTIMIZER 


Used to clear old results table. 



{EditGoto [lidar]EAP:Al } {LET 
guess, le-6} {LET $save_beta, $beta} 

Goto the existing EAP plot (allways 
exists), set the OPTIMIZER guess value 
to a small value and save the current value 
of the backscatter. 

{LET $save_tscope,$tscope_d} {LET 
$save_angle,$align_angle } 

Save exisiting values of the telescope 
diameter and misalignment angle. 

{LET $save_energy, $nrg} {WindowsOff} 

Save existing laser pulse energy and turn 
screen updates off during macro execu- 
tion. 

{EditGoto Cover:Al } {DODIALOG 
EAPoptions,dialog_ans,eap_option } { IF 
dialog_ans=0 } { RETURN } 

Goto a pretty picture! Then display the 
dialog box offering a choice of eap plots. 
Choice is stored in eap_option as a 0 or 1 . 
if user cancels this dialog quit. 

{IF eap_option=l } {setup_varyangle} 

If user chose fixed backscatter and vary- 
ing angle jump to a routine to setup that ; 
calculation. * 

{IF eap_option=0} {setup_varybeta} 

If user chose fixed angle and varying 
backscatter jump to setup routine for that 
option. 

{EditGoto Cover: A 1 } { DODIALOG 
loopprops,dialog_ans,loop_vars } { IF 
dialog_ans=0 } { RETURN } 

Display the input parameter dialog box so 
the user can enter values. If they cancel 
the dialog then quit. 
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B 

c 

{ SELECTBLOCK 

+”EAP:A3 1 ..EAP:”&@CHAR(coI_num 
ber)& @ STRING(row_number+old_nstep 
,0)} (DEL) 

This deletes the results table from the pre- 
vious EAP plot. 

{calc} {LET 

col_number,@CODE(“B”)} {LET 
row_number,32 } 

This resets column and row tracking cells 
used to build results table. 

{FOR loopl_count, 
tscope_min,tscope_max,tscope_step, 
loop2_calc } 

Loop over the telescope values. 

{GraphEdit EAP} {GraphSettings.Type 
“XY,2-D”} 

Finished calculating results so edit EAP 
plot to display results. 

{ Series.Data_Range “XAxisLa- 
belSeries”, 

+”EAP:B3 1 ..EAP:”&@CHAR(col_numb 
er-l)&”31”,l} 

Set up the x-axis range (telescope diame- 
ter). 

{Series.Data_Range “Legend- 
Series”,+”EAP:A32..EAP:A”&@STRIN 
G(row_number+eapvar_nstep,0),l } 

Set up the values for the legends (either 
misalignment angle or beta) 

{FOR 

loopl count, l,eapvar nstep+ 1,1, series ft 

11} 

For each value of misalignment angle or 
beta generate fdl a series with data as a fn. 
of tscope diameter. 

{ Series.Go } { WindowClose } 

Carry out graph changes and close the 
graph edit window. 

{FOR 

loop 1 _count, 1 ,eap var_nstep+ 1,1, series_c 
hng} 

Run a routine to edit each of the data 
series. 

{ SETOB JECTPROPERTY 
“EAP:G$YlAxis.Major_Grid_Style”,”S5 
W 1,0, 0,0”} 

Set up Y-axis grid style. 

{ SETOB JECTPROPERTY 
“EAP:G$Y 1 Axis.Text_Font”,” Anal, 1 0, Y 
es, No, No, No”) 

Set up font for Yaxis labels. 

{ SETOBJECTPROPERTY 

“EAP:G$X 1 Axis.Text_Font”,”Arial, 1 0, Y 

es,No,No,No”} 

Set up font for Xaxis labels. 
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136 series_chng 


B 

c 

{ SETOBJECTPROPERTY 
“EAP:G$Leg- 

end.Text_Font”,”AriaI,10,Yes,No,No,No” 

> 

Set up font for series legends. 

( SETOBJECTPROPERTY 
“EAP:G$YlTitle.Text_Font”,”Arial,12,Y 
es, No, No, No”} 

Set up Yaxis title font. 

{ SETOBJECTPROPERTY 
“EAP:G$XlTitle.Text_Font”,”Arial,12,Y 
es, No, No, No”} 

Set up Xaxis title font. 

(IF eap_option=l } (GraphSettings.Titles 
“EAP plot, constant backscatter 
value”, ’’solved for energy to give 
Pgood=0.5. ’’.Telescope diameter 
(m), Pulse energy (J),””} 

If varying misalignment angle then 
change graph title and axis labels to 
match option chosen. 

(IF eap_option=0} (GraphSettings.Titles 
“EAP plot, constant misalignment 
angle”, ’’solved for energy to give 
Pgood=0.5.”, Telescope diameter 
(m),Pulse energy (J),””} 

If varying backscatter then change graph 
title and axis labels to match option cho- 
sen. 

(LET $tscope_d, $save_tscope}(LET 
$beta,$save_beta} 

Restore telescope diameter and backscat- 
ter to their original values. 

{ LET $nrg,$save_energy } { LET 
$align_angle,$save_angle } 

Restore energy and misalignment angles 
to their original values. 

(LET old_nstep,eapvar_nstep } 

Save value of number of second indepen- 
dent variable steps used. This is needed 
next time the EAP routine is called so that 
the old results table can be cleanly 
deleted. 

(EditGoto EAP:A1} 

Goto the “new” EAP plot and finish. 




(calc } { SETOBJECTPROPERTY 

+”EAP:G$Series[“&@STRING(loopl_c This gives each series a unique marker 
ount,0)&”, l].Marker_Style”,+marker_sty style. 

le» 
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138 series fill 


B 

c 

{ Series.Data_Range loopl_count, 
+”EAP:B”&@STRING(+3 1+loop l_coun 
t,0)&”..EAP:”&@CHAR(col_number- 
1)&@STRING(+31+Ioopl_count,0),l } 

This fills each series with data from the 
results table. 


139 


140 


141 loop2_calc 




{LET $tscope_d,loopl_count} 


{LET 

+”EAP:”& @CHAR(col_number)&”3 1 ”,1 
oopl_count} 


{ IF eap_option= 1 } { FOR 
loop2_count,eapvar_min,eapvar_max,eap 
var_step,eap_calc_angle } 


{IF eap_option=0} {FOR 
loop2_count,0,eapvar_nstep, 1 ,eap_calc_b 
eta} 


{LET col_number,col_number+l } 


{LET row_number,32} 



{ SETOB JECTPROPERTY “loop- 
148 setup_varyangle props:eapvarname.Label_Text”, “Mis- 
alignment Angle”} 


{SETOB JECTPROPERTY “loop- 
props:eapmintext.Label_Text”, “Mini- 
mum (|irad)”} 


{SETOB JECTPROPERTY “loop- 
props :eapmax text. Label_Text”, “Maxi- 
mum (prad)”} 


{SETOB JECTPROPERTY “loop- 
props:eapsteptext.Label_Text”, “No. of 
steps” } 


{LET $EAP:$D$25, +”Backscatter value 
of’} 




This routine is run to set up the loop 
parameters dialog box for varying the 
misalignment angle. 


The labels for the input boxes are changed 
to reflect the required parameters. 



This sets up the “Backscatter value of 
‘value’ /(m-sr). line below the EAP plot. 
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(LET $EAP:$I$15, +”angle (prad).”} 


{ SETOB JECTPROPERTY 
“$EAP:$I$ 14. Alignment”, “Center”} 



158 


{ SETOB JECTPROPERTY 
“$EAP:$I$ 1 5 .Alignment”, “Center” } 

As this is the last line in the subroutine 
execution of the main macro routine 
resumes. 

159 




160 

eap_calc_angle 

{LET $align_angle, +loop2_count} 

This routine is used when misalignment 
angle is the second independent variable. 
It sets up and runs the solver to determine 
the pulse energy required for a given mis- 
alignment angle and telescope diameter. 

161 


{ LET $NRG, Guess } { solver } 

Set energy to a small value. Call the 
solver routine. 

162 


{IF col_number=66 } { LET 

+”EAP:A”&@STRING(row_number,0), 

$align_angle} 

This is one of the routines that creates the 
results table. If this is the first loop over 
the alignment angle write the misalign- 
ment angle in column A. 

163 

' 

{LET 

+”EAP:”&@CHAR(col_number)& @ ST 
RING(row_number,0),$NRG } 

Write the solved value of the energy into 
the results table. 

164 


{LET row_number,row_number+l } 

Increment the row (misalignment angle) 
loop and return to the routine from which 
this one was called. 




{SETOB JECTPROPERTY “loop- 



This routine sets up the labels on the input 


166 setup_varybeta props:eapvarname.Label_Text”, “Back- parameter dialog box for varying back- 


scatter” } 


scatter. 























































178 eap_calc_beta 



{ SETOB JECTPROPERTY “loop- 
props:eapmintext.Label_Text”, “Mini- 
mum /(m-sr)”J 


{SETOB JECTPROPERTY “loop- 
props:eapmaxtext.LabellText”, “Maxi- 
mum /(m-sr)”} 


(SETOB JECTPROPERTY “loop- 
props:eapsteptext.Label_Text”, “No. of 
steps” } 


{LET $EAP:$D$25, +”Misalignment 
angle is”} 


{LET $EAP:$F$25, +$align_angle} 


{LET $EAP:$G$25, +”prad.”} 


{LET $EAP:$I$14, +”backscatter”} 


{LET $EAP:$I$15, +”(/m-sr).”} 


{ SETOB JECTPROPERTY 
“$EAP:$I$ 14. Alignment”, “Center”} 


{ SETOB JECTPROPERTY 
“$EAP:$I$ 15. Alignment”, “Center”} 



This sets up the line below the plot that 
indicates the value of misalignment used. 



This modifies the label on the right of the 
plot to correctly indicate that backscatter 
is the second independent variable. 



Make sure label alignment is correct. 


Return to main routine. 






{LET $beta, +10 A (a+bval*loop2_count)} 

This routine is called when backscatter is 
the second independent variable. This 
selects an intermediate values of back- 
scatter such that values are uniformly dis- 
tributed along a log axis. 

{ LET $NRG, Guess } { solver } 

Set energy to a small value and call the 
solve routine. 

(IF col_number=66} {LET 

+”EAP: A”& @STRING(row_number,0), 

$beta} 

If this is the first loop over the range of 
backscatter values write the backscatter 
value into the first column of the results 
table. 
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B 

c 

{LET 

+”EAP:”&@CHAR(col_number)&@ST 
RING(row_number,0),$NRG } 

Write the solved energy value into the 
results table. 

{LET row_number,row_number+l } 

Increment the row counter and return. 


184 solver 


{ Optimizer.Reset } { Opti- 
mizer.Solution_cell $PBAD } 

This routine solves for the energy 
required. Set OPTIMIZER to solve for 
Pbad. 

{Optimizer.Solution_goal ‘Target 
Value:”} 

Set up as target value for Pbad. 

{Optimizer.Target_Value 0.5 } 

We want Pbad=0.5 (=Pgood). 

{ Optimizer. Variable_cells $NRG} 

We want to vary the energy. 

{Optimizer.Solve} 

Solve! Then return to calling routine. 


Appendix (Id) The expanded MLE lookup table 


mega 

\mi 

0 

m 

ha 

k 

amnia 



5.630 


7.0191 


6.9391 


6.985 


9.848 


10.40 


9.743 


11.21 


12.71 


13.99 


17.70 


17.94 


16.40 


641.0 


52.94 


48. 


4.630 


6.612 


1.0665 


1.1591 


1.35911 


1.479 


1.489 


1.3591 


1.293 


2.633 


2.798 


28.861 


0 


6302 


.0698E-0 


IE 

silon 

Id 

elta 


u 


teiKl 


70.69 


0.874231 


1.198 


1.2045 


1.210 


1.0783 


1.110 


0.92911 


0.91983 


2.4001 0.8838 


4.1771 


SKIESEEjlS] 


.1211 H- 1 1 


3.623 


5.441 


15.71 


14.99 


4.23 


3.839 


4.846 


0.202491 9.47 


0.19923 0.80261 


0.27401 0.5350 


0.653 


1.635 


0.963 


0.6408 


0.61835 


1.015 


1.74 


0.3796 




0.2382 


0.49395 


1.1235 


1.087 


0.3824 


0.2717 


0.4312 


0.05017 


0.21661 


0.23701 


0.2706 


iHH 

■■HKl 


0.1758 


0.189 


1.225 


1.043 


1.130 


1.1531 


1.18431 


3.98 


3.512 


150.0 


11.92 


10.215 


1.2535 


1.573 


11 

0.8853 

4 

5. 


0.5492 

1 

0.18741 


0.905 


0.8760 


0.87923 


0.8835 


4.498 


5.2751 


4.5601 


5.62 


6.231 


4.386 



0.3563 


0.4814 




iHSI 

9.614 

81 

1.1 

la 

2.04531 

0.9001 

■EHEijl 

3.195 

3 

0.1401 


14.003 


7.477 


293.11 


270 



16.6 


0.231 a 


0.77831 


4.09511 


5.63021 


2.568 


2.709 


1.495 


62.54 


680.8 


2.6333 


2.6333 


0.019 


0.049 


0.875 


0.885 


0.779111 


EliEEEEEII 


0.91805 


0.84545 


0.84545 


0.213 


4.284 


3.506 


7.563 


6.5971 


5.980 


10.41 


10.41 


17.45 



1.516310.4648 


1.642610.40951 


1.172 


1.172 


1.146 


msmm 
m mm 


0.1823 


0.1462 


0.1358 


0.1458 


0.146 


0.108 


0.080 


0.1 


0.0761 


0.21171 


0.11711 


0.111 


0.1730 


0.1730 


0.873 

la 

0.4349 

1.21031 

0.3021 


20 


































































































































































































































21 




















Appendix II 


Direct global measurements of tropospheric winds 
employing a simplified coherent laser radar 
using fully scalable technology and technique 


Michael .1. Kavava 

NASA/Marshall Space Flight Center, Code EB54, Huntsville, AL 35812 USA 
michael@lidar.msfc.nasa.gov 

Gary D. Spiers 

Center for Applied Optics, University of Alabama, Huntsville, AL 35899 USA 
garys@photon.msfc.nasa.gov 

Elena S. Lobl 

Hughes STX, 620 Discovery Drive, Huntsville, AL 35806 USA 

Jeff Rothermel 

NASA/Marshall Space Flight Center, Code ES43, Huntsville, AL 35812 USA 
jeff@pcjxr.msfc.nasa.gov 

Vernon W. Keller 

NASA/Marshall Space Flight Center, Code PS02, Huntsville, AL 35812 USA 

ABSTRACT 

Innovative designs of a space-based laser remote sensing “wind machine” are presented. 
These designs seek compatibility with the traditionally conflicting constraints of high scien- 
tific value and low total mission cost. Mission cost is reduced by moving to smaller, lighter, 
more off-the-shelf instrument designs which can be accommodated on smaller launch vehi- 
cles. 


Paper 2214-31, SPIE’s International Symposium on Optical Engineering in Aerospace Sensing, Technical Conference 
2214 on Space Instrumentation and Dual-Use Technologies, session on “Faster, Cheaper, Smaller Space Science Optical Instru- 
ments,” Orlando, FL (6 April 1994). 


1. INTRODUCTION 


Measurements of tropospheric winds from space are highly desired for many NASA, NOAA, DOD, DOE, EPA, and com- 
1 2 

mercial applications. ’ These include global climate change research, improved weather forecast accuracy, optimum aircraft rout- 
ing for passenger safety/comfort and for fuel/time savings, commercial shipping storm avoidance, pollution research and regulation 
enforcement, and military planning. Improved weather forecasting will save lives lost to storms, as well as reduce the number of 
false alarms. Numerous studies have indicated that the optimum measurement approach for winds from space is a pulsed coherent 
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laser radar (CLR) . NASA recently completed dual Phase A and B studies, by GE Astro Space and Lockheed ' , for a full tropo- 
spheric profiling instrument named the Laser Atmospheric Wind Sounder (LAWS). The studies both recommended a pulsed CO 2 
laser CLR system with pulse energy near 20 J, pulse repetition frequency (PRF) near 5 Hz, optical diameter near 1 .6 m, and full- 
time operation (100% orbit duty cycle). However, the projected mission costs were deemed unacceptable for this new era of 
reduced NASA resources, and the projected spacecraft resource requirements exceeded the capabilities of small satellites. NASA’s 
Marshall Space Flight Center (MSFC) has therefore conducted an in-house effort to investigate innovative versions of LAWS 
which are smaller, lighter, less expensive; which consume less power, require less heat removal, fit on small spacecraft and launch 
vehicles; which provide valuable engineering data and space heritage; and which still deliver significant science product consistent 
with NASA’s Mission To Planet Earth (MTPE). The working name for these instruments is Autonomous Earth Orbiting Lidar Util- 
ity Sensor (AEOLUS), after the mythical Greek god of wind. This paper reports on the progress to date of this in-house effort. 

2. WHAT IS AEOLUS? 

The AEOLUS project is a series of point designs of a small, lightweight, low cost, low risk instrument for measuring 
winds from space in regions of high aerosol backscatter. The mission goals are to: 

1) provide valuable scientific information in support of NASA’s MTPE, 

2) demonstrate space operation of key technologies required for a full-scale LAWS, 

3) implement a logical, affordable evolutionary process leading to a full-scale LAWS measurement of global tropospheric 

winds. 

The ground rules we followed are: 

1) to use components that are as close to “commercial off the shelf (COTS)” as possible, 

2) to use technology and a measurement technique that is capable of scaling up to make the full-scale LAWS measure- 

ment, 

3) to design a flexible instrument capable of being accommodated on a variety of platforms in space such as small satel- 

lites, the space shuttle, and space station. 

The instrument design process we followed is shown in Figure 1. 

3. MEASUREMENT TECHNIQUE 

Both the choice of a measurement technique for AEOLUS and subsequent consideration of design trade-offs within the 
chosen technique have confirmed that measurement of wind from space is a complex undertaking that contains numerous and 
sometimes subtle interactions between: 

1) the laser radar instrument parameters, 

2) the available laser radar technology, 

3) the capabilities of the instrument bus or carrier, 

4) the launch vehicle, and 

5) the atmospheric target and its assumed properties. 

When the various requirements of the mission that must be achieved simultaneously are considered from a system per- 
spective, the sensor of choice for this mission is a pulsed CLR. These joint requirements include: 

1) horizontal wind measurement accuracy of about 1 m/s or better, 

2) sufficiently large cross-track measurement swath width for areal coverage, 


3) minimal horizontal wind measurement bias, 

4) horizontal wind measurement resolution of about 100 km or smaller, 

5) best possible vertical resolution, 

6) maximum possible allowed horizontal wind magnitude, 

7) minimal prime power, cooling, and downlink data rate requirements, 

8) minimal instrument mass and volume, 

9) maximum mission duration, 

10) eyesafe fluence exposure, 

1 1) low mission cost. 

Requirement 1 leads to the need for nearly collocated, biperspective line of sight (LOS) wind measurements. Because of 
spatial variability in the horizontal wind field, many LOS measurements are needed in each areal resolution element (e.g., 100 km 
x 100 km square) in order to reduce the contribution to error from undersampling. Since the point on the earth directly below the 
spacecraft moves at 7.3 km/s (350 km orbit height), the 100 km resolution requirement means that the time available for an accu- 
rate LOS horizontal wind measurement in a 100 km square is only about 14 s divided by the number of 100 km squares being 
probed in the cross-track direction. Using our current number of 4 cross-track squares yields about 3.5 s for each horizontal wind 
measurement. (The full-scale LAWS had 10 cross-track squares.) Requirement 2 means the laser beam must be scanned. The scan- 
ning may be continuous or step-stare. Continuous-scanning is preferred by the spacecraft designer who must compensate the angu- 
lar momentum. Since continuous scanning involves changing the beam pointing direction with time, an accurate LOS wind 
measurement must occur over a short time interval due to both the varying spacecraft-earth relative velocity along the LOS, and 
due to the varying probe direction of the laser beam in the wind. Thus the number of laser shots per LOS measurement divided by 
the laser PRF must be small. Pulsed CLR satisfies this by making a measurement with as few laser shots as one. By comparison, 
techniques that employ many laser shots per LOS measurement, such as pulsed noncoherent laser radar (NLR), must use a high 
laser PRF. A higher PRF increases the required laser lifetime (in pulses), increases the on-board data processing rate, and increases 
the data downlink rate unless on-board combining of multiple laser shot data is performed. This is true even if the total transmitted 
pulse energy per LOS wind measurement is the same as CLR. To date, however, NLR measurements have proven less photon effi- 
cient than CLR, and more total energy per LOS measurement must be transmitted. The important comparisons of prime power and 
cooling require knowledge of each laser’s efficiency. Some methods of combining laser shot data on board increase the risk of 
velocity bias, and preclude data processing improvements on the ground. 

4. COMPONENTS OF A SPACE-BASED CLR 

It was stated above that the system designer must account for the interactions of the CLR instrument, the atmospheric or 
earth surface target of interest, the launch vehicle, and the carrier or bus connecting the instrument to the launch vehicle. Each of 
these items has levels of further subdivision. A CLR broadly consists of a laser transmitter, an optics subsystem, a receiver, a con- 
trol computer, a data acquisition subsystem, and a data processing computer. A simplified schematic of a CLR is shown in Figure 2. 
The laser subsystem consists of an optical pulse generating unit and one or two continuous-wave (CW) lasers to perform the func- 
tions of master oscillator and local oscillator (LO). The optics consist of lenses, mirrors, beamsplitters, the beam expanding tele- 
scope, the beam scanner, and, if needed, a lag angle compensation (LAC) element. The receiver comprises an optical detector, and 
optics to combine the received photons from the atmosphere with the LO optical field. This combination must strive to match the 
shape, direction, curvature, and polarization of the LO field with the expected value of signal field. 

5. SPACECRAFT AND LAUNCH VEHICLE ACCOMMODATION 

The payload of the launch vehicle (i.e., spacecraft) consists of the carrier or bus mated to the scientific instrument(s). This 
payload must be accommodated by the launch vehicle in the categories of mass, volume, orbit height, orbit inclination, survivable 
launch vibration frequencies and amplitudes, and survivable launch accelerations. The cost of the launch vehicle and carrier must 
both fit in the mission budget. We are considering Pegasus XL, Conestoga, LLV, and Taurus class launch vehicles. The carrier must 
accommodate the needs of the instrument including electrical power, heat removal (thermal), orbit height, orbit inclination, electri- 
cal power storage in batteries, mass, data downlink rate, mechanical mating, pointing control and stability, etc. 

A sun-synchronous orbit is preferred from the electrical power and heat removal viewpoints. The measurements occur 
over all latitudes and over all the earth after a few days. By contrast an equatorial orbit does not measure the higher latitudes, but 


gives more coverage of low latitudes. It receives the launch advantage of the earth’s rotation, but is less desirable for electrical 
power and heat removal. 

Lower orbit heights allow a greater payload mass, but some of this mass gain is consumed by propellent needed for addi- 
tional reboosts due to increased atmospheric drag. Lower heights also provide greater SNR for fixed nadir angle, but not necessar- 
ily for fixed cross-track coverage. Atmospheric drag also depends on the sunspot cycle, which will peak from 1999 to 2004. The 
designs described later assume a 350 km orbit height with a 30° nadir angle which yields a 400-km wide coverage centered below 
the spacecraft. 

The data rate is proportional to the laser PRF, the number of desired range gates (the height resolution and coverage), and 
the fraction of each orbit spent taking data (orbit duty cycle). The data rate may be lowered by processing data on board before 
transmission. Maximum science is achieved by 100% orbit duty cycle. Downlinking all data allows research into optimum velocity 
estimation algorithms. 

Assuming a suitably located ground station with a large (-5 m) diameter S-band antenna, data can be downlinked directly 
from the spacecraft to the ground station at a peak data rate of ~2 Mbit/s, but only during the time the spacecraft has LOS visibility 
to the receiving station. This available time on any given orbit varies from zero to >5 min. with a gap between downlink opportuni- 
ties sometimes as long as 20 hrs. Thus on-board data storage of several hundred Mbit may be required. A typical ground station 
with a small (~1 m) diameter antenna will handle only a few 10’s of kbit/s. This would require use of several ground stations. At 
greater cost, data could be downlinked through an S-band omnidirectional antenna to the 26-m ground antenna network of the 
Deep Space Network (DSN), or via a geostationary communications satellite such as the Tracking and Data Relay Satellite System 
(TDRSS). 


6. PULSED CLR DESIGN TRADES 

Choosing a design point for a pulsed CLR to measure winds from space leads one to work in an onerous multi-dimen- 
sional space consisting of a dimension for each parameter of the CLR instrument, the atmospheric target, the carrier capabilities, 
and the launch vehicle capabilities. All of these parameters are interconnected and must be dealt with. Some of these interconnec- 

g 

tions have been previously discussed. Varying one parameter to solve one problem often produces one or more new problems. 
Often a mean value description is inadequate and a parameter’s probability distribution function (PDF) must be used (e.g., aerosol 
backscatter). Possible correlations among parameters must also be considered (e.g., aerosol backscatter vs. wind turbulence vs. 
wind shear). Some design choices are continuous (e.g., laser pulse energy, optical diameter), while others are discrete (e.g., velocity 
calibration with or without an earth surface return, conical scanning with a rotating telescope or a rotating wedge or a rotating flat 
mirror, active or passive heat removal). 

An illustrative example of a design trade is the choice of design nadir angle 0, the angle between the nadir direction and 
the LOS laser beam direction, for a laser beam conically scanning about nadir. (Ail our designs assume a constant nadir angle. A 
variable nadir angle would, increase complexity, risk, cost, and mass.) Many aspects of the system design depend on nadir angle. In 
some cases the dependence on nadir angle can be easily pictured, while in other cases the dependence is complex, and should be 
modeled on a computer. Table 1 summarizes some of the simple cases. We assume for the moment a flat earth, and do not distin- 
guish between the nadir angle as the photons leave the CLR in space, and their nadir angle in the atmosphere. This assumption 
becomes more seriously incorrect as the nadir angle increases. Seeking first to maximize signal-to-noise ratio (SNR) (see case 1 in 
Table 1 ), we consider that the signal power falls inversely as the square of the slant range R. The flat earth slant range varies as 
l/cos0. The noise power is proportional to the detection or search bandwidth Bg, which equals (2/A,)Vg, where Vg is the velocity 
search bandwidth. For a maximum design horizontal wind in any direction, we find Bg = (2/X)(2Vjjj^)sin0. The total SNR 

function is therefore (cos 2 0)/sin0. The nonintuitive optimum nadir angle appears to be 0°, or straight down. However, it has been 
shown that SNR is not a good figure of merit (FOM) for wind velocity measurements. Picturing the detected signal in the frequency 
domain, the signal power is gathered in just a few, or even one, frequency bins, while the noise power, assumed white, is spread 
across all frequency bins. A better FOM is the ratio of the signal height to the level of the “grassy” noise, i.e., the ratio of the total 

average signal energy in the observation time to the spectral density level of the noise . This parameter O = SNR x M, where M 

is the number of data samples used to make the frequency estimate^. The number of data samples per estimate is proportional to 


both the sampling frequency of the data recorder and the observation time for one estimate. Holding the desired height resolution of 
the measurement constant means that the observation time goes as l/cos0. Holding Vpjj^ constant, and sampling the data stream 
fast enough to avoid frequency aliasing causes the sampling frequency to go as sin0. This is case 2 in Table 1 and again yields an 
optimum nadir angle of 0°, but with a weaker function of 0. Next we consider that the mission goal is horizontal wind measure- 
ments and not just LOS measurements. The alignment of each LOS measurement with the horizontal wind goes as sin0 and is 
shown in case 3 of Table 1. The optimum nadir angle is now 45°. This combination of LOS velocity estimation performance with 
horizontal alignment is heuristic and not rigorous. Two techniques of effecting a conical scan about nadir are a rotating telescope 
and a rotating wedge. If a wedge is used to deflect the laser beam by 0, the CLR optical diameter is reduced by the factor cos0. 
This effect is added in case 4, where we specifically assume that atmospheric refractive turbulence effects are small, and that there 
is no transmitter/delayed back-propagated local oscillator (BPLO) misalignment angle Ot. The angle Ot occurs between the trans- 
mitted laser pulse direction and the direction of the imaginary BPLO beam as the backscattered photons reenter the CLR. The opti- 
mum nadir angle becomes 35°. We may also heuristically add the goal of large cross-track coverage by the satellite sensor. For a 
flat earth, we multiply by tan0. This is included in case 5. The optimum angle becomes 55°. Case 6 assumes a rotating telescope, 
removes the wedge scanner, and yields that larger nadir angles are always better. We do not use the FOMs in cases 5 and 6 since the 
treatment of swath width is so arbitrary. Our computer simulation for estimating the performance of candidate mission designs uses 
a spherical earth, and also includes the more complex effects of atmospheric extinction, atmospheric refractive turbulence, and mis- 
alignment angle. Using our more sophisticated simulation for cases 4 and 5 produces the smaller optimum angles of 33° and 47°, 
respectively (see Figure 3). It is reasonable that addition of extinction and a spherical earth would favor smaller nadir angles. The 
simulation parameters correspond to our point design 5, which is discussed below. Note that we used a constant misalignment 
angle of 7.3 p,rad, which causes a budgeted 3 dB misalignment loss at our nominal nadir angle of 30°, but varying loss at other 
nadir angles. 

Future refinements to our FOM for just the nadir angle design trade might include the effects of nadir angle on aerosol 
backscatter, land backscatter, ocean backscatter, and the probability of intercepting clouds. Also related are the cost and mass of the 
wedge, and the effects on atmospheric shot spacing of wedge rotation rate and laser PRF. 


Table 3: Selection Of Design Nadir Angle 


Case 

SNR 

(cos 2 0)/sin0 

Number 

Data 

Samples 

tan0 

Horizontal 

Alignment 

sin0 

Wedge 

Scanner 

COS0 

Cross- 

Track 

Width 

tan0 

Total 

Function 

Optimum 

Nadir 

Angle 

(deg.) 

I 

4 





(cos 2 0)/sin0 

U 

2 

4 

4 




COS0 

0 

3 

4 

4 

4 



cos0sin0 

45 

4 

4 

4 

4 

4 


(cos 2 0)sin0 

35 

5 

4 

4 

4 

4 

4 

(sin 2 0)cos0 

55 

6 

4 

4 

4 


4 

sin 2 0 

90 


The misalignment angle a is another important design trade issue, especially when selecting the laser wavelength of the 
8 10 

CLR , or when selecting the optical diameter. Frehlich has examined the effect of misalignment angle a on degradation of the 
CLR heterodyne or mixing efficiency, T|mix; which is important since both SNR and <1> are proportional to it. We employ Fre- 
hlich’s exact results for a circular optical aperture of diameter D, a Gaussian transmitted beam with diameter optimally matched to 
D, a monostatic CLR, far field operation, and negligible refractive turbulence effects. For degradations less than 15 dB, the points 

'y 

approximately follow = 0-42 exp[-(q/2.8) ], where q = (7lDa/X). This is a strong function of the ratio q, which includes 
optical diameter, misalignment angle, and optical wavelength. If a wedge scanner is used, the diameter D refers to the smaller value 



exiting the wedge. For a fixed budgeted loss in SNR due to misalignment, the ratio q must be held constant. Smaller values of A. 
require proportionally smaller values of allowed misalignment. This is shown in Figure 4 for two wavelengths and the parameters 
of our point design 5. When Ct = 0 and effects of refractive turbulence are negligible, the common intuition that larger values of D 

yield better performance is correct since SNR is proportional to receiver area = 7tD/4. However, nonzero values of CL cause a 
reduction in this quadratic gain of performance with increasing D, and may even cause an absolute reduction in performance. This 
is shown in Figure 5 using the parameters of our point design 5. Even without the misalignment effect, thus having the promise of 
quadratic gain in performance, the optical diameter would be limited due to the penalties of mass, volume, cost, need for high opti- 
cal quality, and the need to conically scan the beam. Misalignment even further reduces the optimum design diameter. We are 
attempting to merge the behavior of the CLR with information about realistic on-orbit misalignment angles. The actual misalign- 
ment angle will consist of contributions from prelaunch assembly and alignment of the optical subsystem, further misalignment 
from launch stress and orbit life, laser shot pointing jitter, and spacecraft pointing jitter during the round trip time of the transmitted 
photons (- 3 ms). Techniques to eliminate the first three contributors are possible, but they complicate the instrument design. 

7. INSTRUMENT DESIGNS 

The NASA/MSFC AEOLUS team has nearly completed five instrument point designs. We are carrying two candidate 
laser technologies, the CO 2 laser at 9.1 1 fim wavelength, and the Tm,Ho:YLF diode-pumped solid-state laser at 2.06 p.m. From 
our nearly COTS ground rule, and our desire for smallsat capability, we have permitted maximum laser pulse energies of 400 and 
200 mJ, respectively, and a maximum optical diameter of 50 cm. Each design is examined from the mechanical, electrical, thermal, 
optical, laser, CLR wind measurement performance, and spacecraft accommodation perspectives. An example of CLR perfor- 
mance prediction for our point design number 5 is shown in Figure 6. Five curves are plotted against the atmospheric aerosol back- 
scatter coefficient. The reflectances of land and water may also be converted into these backscatter units. This description permits 
scientists to determine the applicability of each CLR to measure winds of various atmospheric targets such as clear air, clouds, dust, 
boundary layer air, jet streams, etc. From top to bottom, the five curves are SNRj^j, the unrealistic limiting case of matched filter 

SNR, SNRg, the more realistic SNR calculated from the noise admitted into the velocity search bandwidth, <D, as discussed earlier, 

o 

Pg, the probability that a LOS wind estimate is “good”, and CTy, the standard deviation or spread of the “good” wind estimates . 
The “good” wind estimates are clustered around the true value of the wind speed, while the “bad” estimates are uniformly distrib- 
uted over a region of width Vg. All the values of (Jy are sufficiently small to provide excellent performance. Even as Pg falls to 
20% for low values of backscatter, CTy remains smaller than 0.4 m/s. Therefore, the key performance criterion to monitor is Pg. 

Since the backscatter value that yields P g = 0.5 occurs near P(2.06 |im) = 2.5 x 10’^ m'^sr’*, we refer to that number as the aerosol 
backscatter sensitivity of design 5. As can be seen, lower values of P(2.06 jtm) will produce fewer but usable velocity estimates. 
Note that for P(2.06 |im) = 2.5 x 10"^ m**sr'*, SNRj^j - 3 dB, SNRg - -15 dB, O ~ 6, and ay - 0.3 m/s. 

All five of our point designs are contrasted in Table 2. The top row shows the backscatter sensitivity of each design. Care 
must be taken when comparing these values at different optical wavelengths, since backscatter varies with wavelength. Common to 
all 5 designs are an SNR margin of 3 dB, a budgeted misalignment loss of 3 dB, a sun-synchronous orbit height of 350 km, a 30° 
nadir angle, a target altitude of 300 m, transmit and receive optics efficiencies of 0.9 each, and no polarization mismatch loss. Fig- 

7 

ure 7 contrasts the sizes of designs 1-5 with the Lockheed LAWS Phase B design . Design 5 has a pulse energy-receiver aperture 
product 15 dB below designs 3-4 in order to be accommodated on a Pegasus rocket It achieves considerable volume and mass sav- 
ings. The orbit average electrical power needed by design 5 is higher than 3-4 due to the increase of orbit duty cycle to 30% (a fac- 
tor of 6), and an increase in laser PRF to 50 Hz (a factor of 5). Designs 2,4, and 5 use a rotating wedge to create a conical scan 
about the nadir direction. Designs 1 and 3 have two fixed pointing directions, fore and aft, to allow biperspective wind measure- 
ments along a line parallel to the ground track, but offset by 144 km. The two views of a single point in the atmosphere would be 
accomplished by switching the CLR between two 50 cm telescopes, and would occur about 40 s apart. This configuration has less 
science value than the conical scan. Comparing the 9. 1 1 J4m designs 1-2 with the 2.06 |im designs 3-4, the cost of scanning in mass 
and power, as well as its lower volume can be seen. 

No attempt has been made to have equal science value between the point designs at 2.06 and 9.1 1 Jim. A study is under- 
way to quantify the ratio of backscatter values between the two wavelengths for various candidate atmospheric targets within the 


sensitivity range of our different point designs. Note also that holding the misalignment loss fixed at 3 dB causes the misalignment 
angle specification to become stricter at 2.06 fj.m. Our designs are not yet mature enough to quantify the effect this stricter specifi- 
cation might have on instrument mass, volume, complexity, and cost; or on the suitability of candidate spacecraft. 

8. CONCLUSIONS 

Global tropospheric wind measurements are highly desired and will provide many benefits. Numerous studies have 
selected coherent laser radar as the optimum technique. Atmospheric winds have been successfully measured with coherent laser 
radar since 1967. A full-scale mission which measures the lowest levels of aerosol backscatter is not practical in today’s economic 
climate. However, significant scientific benefits are possible with a smaller and more affordable instrument. Significant benefit 
could result in only a few years by starting a faster, smaller, cheaper mission now. NASA/MSFC is ready to design and perform 
such a mission. The mission will provide exciting science and will also provide valuable information for a future full-scale effort. 
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DESIGN 

PARAMETER 

NO. 1 

NO. 2 

NO. 3 

NO. 4 

NO. 5A 
THIN CLOUDS 

NO. 5B 
THICK CLOUDS 

PERFORMANCE 







BACKSCATTER(/M-SR)«50% 

8.2E-09 

1.1E-08 

5.7E-08 

7.6E-08 

2.5E-06 

3.9E-05 

SNR MARGIN, dB 

3 

3 

3 

3 

3 

3 

VERTICAL RESOLUTION, M 

1000 

1000 

1000 

1000 

1000 

64 

MAXIMUM HORIZONTAL WIND. M/S 

30 

30 

30 

30 

50 

50 

INSTRUMENT 







LASER WAVELENGTH, pM 

9.11 

9.11 

2.06 

2.06 

2.06 

2.06 

SCAN TYPE 

dual-look 

scan wedge 

dual-look 

scan wedge 

scan wedge 

scan wedge 

PULSE ENERGY, mJ 

400 

400 

200 

20 0 

25 

25 

PRF.HZ 

20 

20 

10 

10 

50 

50 

TELESCOPE DIAMETER, M 

0.5 

0.5 

0.5 

0.5 

0.25 j 

0.25 

MISAUGNMENT LOSS, dB 

3 

3 

3 

3 

3 

3 

MISALIGNMENT ANGLE. pRAD 

13.9 

16.1 

3.15 

3.64 

7.28 

7.28 

SIZE, M 

1.22x1.24x1.33 

1.22(D)x1.35 

1.22X1.24X1.33 

1.22(D)x1.35 

0.73(D)x0.99 

0.73(D)x0.99 

Inches 

48.0x48.6x52.4 

48.0(D)x53.3 

48.0x48.8x52.4 

48.0(D)x53.3 

28.8(D)x38.9 

28.8(D)x38.9 

VOLUME, M3 

2 

1.6 

2 

1.6 

0.4 

0.4 

ft3 

71 

55.8 

71 

55.8 

14.7 

14.7 

MASS, Kg 

231 

266 

191 

226 

125 

125 

lbs 

508 

585 

420 

497 

275 

275 

POWER, W (STBY/WARM-UP/OPER) 

160/332/437 

168/378/476 

75/182/355 

83/226/391 

83/221/331 

83/221/331 

ORBiT DUTY CYCLE, % 

5 

5 

5 

5 

30 

30 

ORBIT AVERAGE POWER. W 

203 

219 

107 

122 

180 

180 


Table 2 AEOLUS designs intercomparison. 






















Figure of Merit (A.U.) 



Figure 1 AEOLUS instrument design process. 



Figure 3 Two selected figures of merit as a function of nadir angle (see text). 
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Figure 2 Optical Schematic for AEOLUS point design no. 5. 
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Figure 4 Heterodyne efficiency vs. misalignment angle for a 25 cm telescope, wedge scanner and 30 degree nadir angle. 
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Figure 5 The variation of matched filter SNR with telescope diameter and misalignment angle for a wavelength of 2pm. 
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Figure 6 A representative performance analysis for Design No. 5 with a 1 km vertical integration range. 



LAWS Phase B 


AEOLUS AEOLUS AEOLUS 

Designs 1 & 3 Designs 2 & 4 Design 5 


Figure 7 Comparison of AEOLUS designs and one of the LAWS Phase B designs. 





